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Abstract 

In this paper matrix orthogonal polynomials in the real line are described in terms of a 
Riemann-Hilbert problem. This approach provides an easy derivation of discrete equations 
for the corresponding matrix recursion coefficients. The discrete equation is explicitly derived 
in the matrix Freud case, associated with matrix quartic potentials. It is shown that, when 
the initial condition and the measure are simultaneously triangularizable, this matrix discrete 
equation possesses the singularity confinement property, independently if the solution under 
consideration is given by recursion coefficients to quartic Freud matrix orthogonal polynomials 
or not. 

1 Introduction 

The study of singularities of the solutions of nonlinear ordinary differential equations and, in partic- 
ular, the quest of equations whose solutions are free of movable critical points, the so called Painleve 
property, lead, more than 110 years ago, to the Painleve transcendents, see pQ (and [2] for a recent 
account of the state of the art in this subject). The Painleve equations are relevant in a diversity of 
fields, not only in Mathematics but also, for example in Theoretical Physics and in particular in 2D 
Quantum Gravity an Topological Field Theory, see for example [2]. 

A discrete version of the Painleve property, the singularity confinement property, was introduced 
for the first time by Grammaticos, Ramani and Papageorgiou in 1991 when they studied some 
discrete equations, including the dPI equation (discrete version of the first Painleve equation), see 
also the contribution of these authors to [2j. For this equation they realized that if eventually a 
singularity could appear at some specific value of the discrete independent variable it would disappear 
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after performing few steps or iterations in the equation. This property, as mentioned previously, is 
considered by these authors as the equivalent of the Painleve property [I] for discrete equations. 
Hietarinta and Grammaticos also derived some discrete versions of the other five Painleve equations 
[U E]. See also the interesting papers [6] and [7j. 

Freud orthogonal polynomials in the real line [8] are associated to the weight 

w p (x) = |x| p e" |x|m , p > -1, m > 0. 

Interestingly, for m — 2,4, 6 it has been shown [U] that from the recursion relation 

xPn = a n+1 p n+1 (x) + b n p n {x) + a n p„_i(a;), 

the orthogonality of the polynomials leads to a recursion relation satisfied by the recursion coefficients 
a n . In particular, for m = 4 Van Assche obtains for a n the discrete Painleve I equations, and therefore 
its singularities are confined. For related results see also [10]. For a modern and comprehensive 
account of this subject see the survey [TT] . 

In 1992 it was found [12] that the solution of a 2 x 2 Riemann-Hilbert problem can be expressed 
in terms of orthogonal polynomials in the real line and its Cauchy transforms. Later on this property 
has been used in the study of certain properties of asymptotic analysis of orthogonal polynomials 
and extended to other contexts, for example for the multiple orthogonal polynomials of mixed type 

m. 

Orthogonal polynomials with matrix coefficients on the real line have been considered in detail 
first by Krein [T4"l [TB] in 1949, and then were studied sporadically until the last decade of the XX 
century. These are some papers of this subject: Berezanskii (1968) [IE], Geronimo [H] (1982), and 
Aptekarev and Nikishin [18] (1984). In the last paper they solved the scattering problem for a 
kind of discrete Sturm-Liouville operators that are equivalent to the recurrence equation for scalar 
orthogonal polynomials. They found that polynomials that satisfy a recurrence relation of the form 



xP k (x) = A k P k+l {x) + B k P k {x) + A? k _ x P k .. x {x), k = 0, 1 



are orthogonal with respect to a positive definite measure. This is a matricial version of Favard's 
theorem for scalar orthogonal polynomials. Then, in the 1990's and the 2000's some authors found 
that matrix orthogonal polynomials (MOP) satisfy in certain cases some properties that satisfy 
scalar valued orthogonal polynomials; for example, Laguerre, Hermite and Jacobi polynomials, i.e., 
the scalar-type Rodrigues' formula [191 12Q1 [21] and a second order differential equation [221 [231 [21] . 

Later on, it has been proven [22] that operators of the form D=d 2 F 2 (t)+d 1 Fi(t)+d° F have as 
eigenfunctions different infinite families of MOP's. Moreover, in [21] a new family of MOP's satisfying 
second order differential equations whose coefficients do not behave asymptotically as the identity 
matrix was found. See also [26J . 

The aim of this paper is to explore the singularity confinement property in the realm of matrix 
orthogonal polynomials. For that aim following [12] we formulate the matrix Riemann-Hilbert prob- 
lem associated with the MOP's. From the Riemann-Hilbert problem it follows not only the recursion 



relations but also, for a type of matrix Freud weight with m = 4, a nonlinear recursion relation (58) 
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for the matrix recursion coefficients, that might be considered a matrix version -non Abelian- of the 
discrete Painleve I. Finally, we prove that this matrix equation possesses the singularity confinement 
property, and that after a maximum of 4 steps the singularity disappears. This happens when the 
quartic potential V and the initial recursion coefficient are simultaneously triangularizable. It is im- 
portant to notice that the recursion coefficients for the matrix orthogonal Freud polynomials provide 
solutions to (58) and therefore the singularities are confined. A relevant fact for this solution is that 
the collection of all recursion coefficients is an Abelian set of matrices. However, not all solutions of 
(58) define a commutative set; nevertheless, the singularity confinement still holds. In this respect 
we must stress that our singularity confinement proof do not rely in matrix orthogonal polynomials 
theory but only on the analysis of the discrete equation. This special feature is not present in the 
scalar case previously studied elsewhere. 

The layout of this paper is as follows. In section 2 the Riemann-Hilbert problem for matrix 
orthogonal polynomials is derived and some of its consequences studied. In §3 a discrete matrix 
equation, for which the recursion coefficients of the Freud MOP's are solutions, is derived and it is 
also proven that its singularities are confined. Therefore, it might be considered as a matrix discrete 
Painleve I equation. 



2 Riemann-Hilbert problems and matrix orthogonal poly- 
nomials in the real line 

2.1 Preliminaries on monic matrix orthogonal polynomials in the real 
line 

A family of matrix orthogonal polynomials (MOP's) in the real line [11] is associated with a matrix- 
valued measure p on M; i. e., an assignment of a positive semi-definite N x N Hermitian matrix p(X) 
to every Borel set X C IR which is countably additive. However, in this paper we constraint ourself 
to the following case: given an N x N Hermitian matrix V(x) = (Vy(a;)), we choose dp = p(x)dx, 
being dx the Lebesgue measure in K, and with the weight function specified by p = exp(— V{x)) 
(thus p is a positive semi-definite Hermitian matrix). Moreover, we will consider only even functions 
in x, V(x) = V(—x); in this situation the finiteness of the measure dp is achieved for any set of 
polynomials V i; j(x) in x 2 . Associated with this measure we have a unique family {P n (x)}^ =0 of 
monic matrix orthogonal polynomials 

P n {z) = l N z n + t! 1 ^"" 1 + • • • + i { n ] e C NxN , 

such that 

/ P n (x)xi p(x)dx = 0, j — 0, . . . , n — 1. (1) 

Jr 

Here Ijv denotes the identity matrix in C NxN . 
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In terms of the moments of the measure dfi, 

x j p(x)dx G C NxN , 
we define the truncated moment matrix 



ni 



'3 •" 



J =0,1, 



m W ■= (nnj) G C 



nNxnN 



with mij=mi + j and < i,j < n— 1. Invertibility of m^ n \ i.e. detm*™) ^ 0, is equivalent to the 
existence of a unique family of monic matrix orthogonal polynomials. In fact, we can write ([!]) as 



/ m •• 

ym n _i • • • m 2 „- 2 / 



\ 








/ 


Uv 




\-m 2n -i ) 



(2) 



and hence uniqueness is equivalent to detm^ ^ 0. From the uniqueness and evenness we deduce 
that 



P n (z) = I N z n + ^z n - 2 + ^z n ~ 4 + • • • + 7 i n) , 



(3) 



where 7„ = if n is odd. 

The Cauchy transform of P n (z) is defined by 

Qn{z) :-- 



p(x)ax, 



2tt i J R x — z 

which is analytic for z G C\R. Recalling ^ = \ YTjZl fj + \ j/te an d ^ we get 

1 1 



(4) 



and consequently 



Q n (z) = c; 1 ^- 1 + 0(^ 1 - 2 ) 



z — >■ oo, 



where we have introduced the coefficients 



1 



2tt i 



P n (x)p(a;)a; n dx 



(5) 



(6) 



(7) 



relevant in the sequel of the paper. 
Proposition 1. We have that c n satisfies 



det(m (n) ) 
detc B = -27ri det(m(n+1)) . 
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Proof. To prove it just define m := (m n , ...,m 2n -_i), consider the identity 



m (n) 1 0^ m(n+ i) _ (l nN ml"' l m 



Iff J V m< m 2n 

and apply the Gauss elimination method to get 

det(m( n+1 )) , . , , s-i . . n 

- —, — = det(m 2n - rrrmS > m) 0; 
det(mW) v ; r ' 

from ([2]) we conclude m2„ — m'm'"' 1 m = J R P n (x)x n p(x)dx. □ 

The evenness of V leads to Q n (-2 ; ) = (— l) n+1 Qn(— 2), so that 

00 

Q n (z) = c;^-"- 1 + J2 a^- l) z- n - 2j+ \ z -> 00. (9) 

3=2 

In particular, 

Q (z) = c^ 1 z~ 1 + c[ 1 z- 3 + 0(z- 5 ), z^oo. (10) 

Finally, if we assume that are Holder continuous we get the Plemelj formulae 

(Q n {z) + -Q n {z)\ =P n (x)p(x), (11) 

with Q n (x)+ = Q n (z)\z=x+io+ and Q n (x)~ = Q n (z)\ z=x+i0 - . 

2.2 Riemann— Hilbert problem 

Definition 1. The Riemann-Hilbert problem to consider here is the finding of a 2N x 2N matrix 
function Y n (z) E C 2Nx2N such that 

1. Y n (z) is analytic in z G C\R. 

2. Asymptotically behaves as 

_!,v (I N Z n 



^nW = (I 2 V + 0(z- i ))( ^ ^_ n ), *->00. (12) 



5. On K we /jane t/je jump 



y n (x) + = y n (x)_ . (13) 



1 



N 



An easy extension of the connection among orthogonal polynomials in the real line with a par- 
ticular Riemann-Hilbert problem discovered in [12] can be proven in this matrix context. 

Proposition 2. The unique solution to the Riemann-Hilbert problem specified in Definition^is given 
in terms of monic matrix orthogonal polynomials with respect to p(x)dx and its Cauchy transforms: 



. Pn(z) Qn(z) 

Y n (z) =[ ) , n > 1. (14) 

—l{z) C n _iQn-l \ z ) 



/ F n {z) Qn(Z) \ 

Proof. In the first place let us show that ( J does satisfy the three conditions 

requested by Definition [TJ 



1. As the matrix orthogonal polynomials P n are analytic in C and its Cauchy transforms are 
analytic in C\R, the proposed solution is analytic in C\R. 

2. Replacing the asymptotics of the matrix orthogonal polynomials and its Cauchy transforms we 

apt ( Pni - z) Q " {z) \ -V ( zn +0(z n - 1 ) Oiz-"- 1 ) \_ (J[ rX^w f InZ " \ xtrhem 

§ et [c^iP^z) Cn-iOn-iW ) { O(z^) z-+0(z— 1) ) ~ W + UK* ))[ o I N z-) When 

z — >■ oo. 

^ Prnmllnlwpcrptf Pn(x+iO) Q„(x+iO) \_/ P„(z-iO) Q„(x-iO) \ _ / P n (x)p(x) 

o. num <|±jj/ we get ^ Cn _ 1 p ji _ l(:c+i o) cn-iQ«-i(a:+iO) ^ I Cn _iP n _i(x-iO) Cn _iQ„_i( :c -iO) ^ k c n - 1 P n - 1 (x) P (x) J ■ 

Then, a solution to the RH problem is Y n — ( c fp^ l(2 ) c But ^ ne solution is unique, as 

we will show now. Given any solution Y n , its determinant detY n (z) is analytic in C\R and satisfies 

detr n (a;) + = det^F„(a;)_ ^ P ^ \ = detr„(x)_det ^ p ^ 
= detF n (x)_. 

Hence, detY n (z) is analytic in C. Moreover, Definition [l] implies that 

detF n (^) = 1 + 0(^ 1 ), z -> oo, 

and Liouville theorem ensures that 

detY n (z) = 1, \/z E C. (15) 



From (15) we conclude that Y~ x is analytic in C\R. Given two solutions Y n and Y n of the RH 
problem we consider the matrix YnY^ 1 , and observe that from property 3 of Definition [l[ we have 
(Y n Y~ 1 ) + = (Y^Y^ _1 )_, and consequently YnY" 1 is analytic in C. From Definition [I] we get YnY^ 1 — >■ 
I2JV as z — ?• 00, and Liouville theorem implies that YrX~ x = I2N] i-e., Y n = Y n and the solution is 
unique. □ 
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Definition 2. Given the matrix Y n we define 

SJz) := Yjz) 



l N z~ n 
I N z n 



SJz) 



Proposition 3. 1. The matrix S n has unit determinant: 

det SJz) = 1. 

2. It has the special form 

AJz 2 ) z-'Bjz 2 )" 
^CJz 2 ) Djz 2 ) 

3. The coefficients of S n admit the asymptotic expansions 
AJz 2 ) =I N + S$ lZ - 2 + 0(z- 4 ), Bjz 2 ) = Sg 2 + S%z~ 2 + 0(z" 
CJz 2 ) = + S%z~ 2 + 0{z-% Djz 2 ) =I N + S%z- 2 + 0{z~% 

for z — > oo. 



Proof. 1. Is a consequence of (15) and (16) 



2. It follows from the parity of P n and Q n . 



3. (12) implies the following asymptotic behaviour 

SJz) = I 2N + S^z- 1 + 0(z- 2 ), 



and (18) gives 



<j(2i) 



s 



(2i) 
n,ll 



s. 





(20 I > 

n,22 . 



o(2i-l) 



and the result follows. 



z — > oo, 



5 



(2i-l) 
n,21 



o(2i-l) X 
°n,12 



Observe that from (18) we get 

o—i( \ _ ( AJz 2 ) z 1 Bjz 2 

n [z) ~ Kz-'cjz 2 ) bjz 2 ) 

with the asymptotic expansions for z — > oo 

AJz 2 ) =I N + (S%S { ^ 21 - S%)z~ 2 + 0{z'% 

Bjz ) = —S^2 — (5^12 ~ ^n,ll^n,12 + S n ^ 12 (S n^lS nj i2 ~ ^i.L)) Z + @{ z ) 

CJz ) = —5^2! + ( — 5^21 + S^iS^u + (5^22 — S n piS^ 12)^,21) z + 0{z 

bjz 2 ) = i N + (sSis® - s%)z~ 2 + o(z-*). 



2.2.1 Recursion relations 

We now introduce the necessary elements, within the Riemann-Hilbert problem approach, to derive 
the recursion relations and properties of the recursion coefficients in the context of matrix orthogonal 
polynomials. 

Definition 3. We introduce the matrix 

V W \c n ^P n ^{z)p{z) c n ^Q n ^(z)J 
Proposition 4. 1. Z n (z) is analytic on 
2. for z — > oo it holds that 



Z n (z) = (I 2N + 0(z- 1 ))^ 



z n p(z) 
z~ n I N 



i 



3. over K. it is satisfied 



Definition 4. We introduce 



(23) 



Z„(,) + /. ( \ \ N n ). (21) 



M B (s) := ^W(z), (25) 
i? n (z) := ^(z)^" 1 ^) = Y n+l {z)Y n -\z). (26) 

We can easily show that 
Proposition 5. TTie matrices M n and R n satisfy 

M n+1 (z)R n (z) = ^Rn(z) + Rn{z)M n {z). (27) 
Proof. It follows from the compatibility condition 

where TF„ := F n+1 . □ 
We can also show that 
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Proposition 6. For the functions R n (z) and M n (z) we have the expressions 



Rn(z) 



MJz) 



zI N 

°n+l,21 



-s, 



(1) 

n,12 



A n {z^p-\z)A n {z*) A n {z*)z-^p-\z)B n {z 2 ) 



^C n {z*) A -^p-\z)A n {z 2 ) z -2 Cn (z 2 )^fp-\z)B n (z*), 
where [•]+ denotes the part in positive powers of z. 

Proof. The expression for R n is a consequence of the following reasoning: 

1. In the first place notice that R n (z) is analytic for z G C\R. 

2. Moreover, denoting 

R n+ (x) := Y n+ i + {x)(Y n+ (x)) , 



(28) 
(29) 



R n _(x) :=Y n+1 _(x)(Y n _(x)) 



-i 



(30) 
(31) 



and substituting (13) in (30) we get R n+ (x) = R n _(x) and therefore R n (z) is analytic in C. 



3. Finally, if we substitute (16) in (26) we deduce that 



R n (z) = Y n+1 (z)Y n - L (z) 

'zl N 
z- 1 ! 



S n +i( z ) 



N 



s-Hz) 







zl N 0\ „(i) (I N 

o o r n+l I o z-% 



In 
>i N / v 

and the analyticity of R n in C leads to the desired result. 
For the expression for M n we have the argumentation 

1. M n (z) is analytic for z e C\R. 

2. Given 

M n+ (x) := ^^(Z n+ (x))-\ 
M n _(x) := ^^(Z^ix))- 1 . 



SP + 0(z- x ), 



Z — > OO, 



dz 



(32) 
(33) 



Substituting (24) in (32) we get 



M n+ {x) = M n _(x), 
and therefore M n (z) is analytic over C. 
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3. From (16) and (22) we see that Z n (z) is 



so that 



where 



ZJz) = Sjz) 



z n p(z) 
z~ n l 



N 



dZ n (z) z _ 1(z) = ^Ms^z)- 1 + S n (z)K n (z)S-\z), 



dz 



dz 



-nz- 1 ! 



Finally, as M n (z) is analytic over C, (35) leads to 



M n (z) = ^Z-'W 



Observe that the diagonal terms of M n are odd functions of z while the off diagonal are 
functions of z. Now we give a parametrization of the first coefficients of S in terms of c n . 



Proposition 7. The following formulae hold true 



c(i) _ r -i 



o(2) 
°n,ll 



7(3) 



^n,21 — r 



i=l 

n-1 

n— 1 



-c 1 



c(2) _ f . . 

°n,22 ~ i-1 i ' 

i=l 

n+1 

c(3) r- 1 
i=l 



-1 



i=l 



Proof. Equating the expressions for Y n (z) provided by (14) and (16) we get 



Yjz) 



Pn(z) Q n {z) 
Cn-lP n -l\Z) C n -iQ„-i(z) 



z — > oo. 
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Expanding the r.h.s. we get 



J n,21 
y(2) 



c(l) _ r -l 



qw _ n 



°1,21 ~~ °2,21 ~~ u ' °n,21 ~~ c n-l^ n -l 



11' 



c( 3 ) _ 

°n,12 ~~ c n °ra+l,22) 



-lo( 2 ) 



where we have used that 



q(2) _ _ _-i 

°1,22 — c O c 'l j 



which can be proved from (10). Introducing (37) into (28) we get 

RJz) = 













and (16) and (39) lead to 



S n +i{ z ) 



zI N 










SJz) 



'z- l I N 
zl 



so that 



q(2) _ q(2) 
°n+l,ll J i 



-1, 



ri(3) — i ri, , 

*->n,12 ~ C n ^n,22 ~ C n+li 



11 ~~ C n C n-li 
lo(2) _ _-l 



where we have used (37). From (37) and (41) we get 



rf(2) r,(2) _ -1 

' 5 n+l,22 D n,22 — C n- C n+1- 



Summing up in n in (40) and (42) we deduce 

n-l 

E,q(2) _ Q (2) 

Wi+1,11 °i,l] 



n-l 

: ~ 5Z C * 1C * -1 ' 
t=l i=l 

n— 1 n— 1 

V(<? (2) -9^1- Vr-r- 1 

/ yWi+1,22 °i,22y | ~~ / ,, c t c i+l' 



t=l 



i=l 



which leads to 



q(2) 

°n,ll 



J2 C * %Ci - 1 + °n 



C-n—li 



i=l 



s. 



(2) 
n,22 



E 

i=l 



Ci-XCa 
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where we have used (37) and (38). Finally (37), (43) and (44) give 



n-1 



S 



(3) 
n,2l 



C n -2) 



(45) 



i=l 



valid for n > 2, and 



n+l 

c(3) _„-l^„ „- 



i=l 



Notice that (37) gives 



c(3) _ n 



(46) 

□ 

(47) 



Proposition 8. Matrix orthogonal polynomials P n (and its Cauchy transforms Q n ) are subject to 
the following recursion relations 



P n +i{z) = zP n (z) - ]^[5 n P n ^i{z) 



with the recursion coefficients (3 n given by 



n> 1, 



Po ■= 0. 



Proof. Observe that ( 26 ) can be written as 

Y n+1 (z) = R n (z)Y n (z) 



Then, if we replace (14) and (39) into (50) we get the result. 



(48) 

(49) 

(50) 

□ 



We now show some commutative properties of the polynomials and the recursion coefficients. 
Proposition 9. Let f(z) : C ->• C NxN such that [V(x), f(z)\ = V(x, zjeRxC, then 



[c n ,f{z)] = \p n ,f(z)]=0, 
[P n (z')J(z)]=0, 



n > 0, 
n > 0, 



Wz e C, 
Vz, z' e C. 



Proof. Let us suppose that for a given m > we have 

[P m (x),f(z)] = [P m -i(x),f(z)] = 0. 
Then, recalling ^ these expressions give 

[CmJ(z)] = [Cmr-l,f(z)] = 0, 



(51) 



(52) 
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respectively. Therefore, using the recursion relations (48) and (|49|) we obtain 



This means that 



[P m+ i(x)J(z)] = x[P m (x),f(z)] - [c m 1 c m _iP m _i(x),/(2;)] = 0. 



[Cm+1,/0)] = 0. 



(53) 



Hypothesis (51) holds for m=l, consequently [c n ,f(z)] = for n = 0,1,... and (49) implies 
[f3 n ,f(z)] = 0. Finally, as the coefficients of the matrix orthogonal P n {z) are polynomials in the 
/3's we conclude that [P n (z'), f(z)] = for all z, z' G C. □ 



Corollary 1. Suppose that [V(x), V(z)] = for all x 6R and z£C, then 

[P n (z),P m (z')]=0, Vn,m>0, z,z'eC, 

\pm Cm~\ 0, 
\P n ,Pm]=0. 



(54) 
(55) 
(56) 



Proof. Applying Proposition [9] to / = V we deduce that [P n (z'),V(z)} = 0, so that it allows to use 
again Propositi on [9~| but now with f = P n and get the stated result. From (|7| and (54) we deduce 
(55) and using (149 \ we get ([561. □ 



3 A discrete matrix equation, related to Freud matrix or- 
thogonal polynomials, with singularity confinement 

We will consider the particular case when 

V(z) =az 2 + I N z 4 , a = a t . (57) 

Observe that [V(z), V'(z')] = for any pair of complex numbers z, z' . Hence, in this case the corre- 
sponding set of matrix orthogonal polynomials {P n }^L , that we refer as matrix Freud polynomials, 
is an Abelian set. Moreover, we have 



[c n , c m ] = [/3 n , p m ] = [c n , a] = [p n , a] = 0, Vn, m = 0, 1, 

In this situation we have 



Theorem 1. The recursion coefficients /3 n (49) for the Freud matrix orthogonal polynomials deter- 
mined by (57) satisfy 

(3 n+1 = n/3' 1 - f3 n „i - p n -a, n = 1, 2, . . . (58) 

with /3q := 0. 
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Proof. We compute now the matrix M n , for which we have 



MJz) 



and is clear that 



with 



-A n {z 2 ){2az + Az 3 I N )A n {z 2 ) -A n {z 2 ){2a + Az 2 I N )B n {z 2 ) 
-C n (z 2 )(2a + Az 2 l N )A n (z 2 ) -C n (z 2 )(2a Z - 1 + 4zI N )B n (z 2 ] 



M n (z) = M^z 3 + M^z 2 + M^z + Jiff, 



(59) 



(60) 



M (3) = /-4Iiv 0i ( 2) 







4 ^!-2| 



(1) qW 

n,12°n,21 

4S^lS^ } ' ' 





-2 9 (1) o/-49 {3) + 49 {1) 9 (2) -49 (1) 9 {1) 9 (1) 

^" D n,21 LX ^°n,21 ' ^°n,21°n,n ^°n,21 n,12°n,21 



2a ^nA2 + 45^12 + 4^1,12(^,21 ^n,12 ~~ ^,22) 





Replacing (37)- (46) into (60) we get 



()/' ~ \-4c n ^ h M « -{ 4c n _ 1 c; 1 '' 



M{ 0) 
Mi°) 







—4cqC x Cq — 2c a 


— 4c n _2 — 4c ra _iC n 1 c n _i — 2c n _i« 



4c 2 1 + 4q 1 c c 1 1 + 2ac x 1 


4 c n+i "I" 4c n c n _ lCn 1 4"2«C n 1 







(62) 

n > 2. 
(63) 



The compatibility condition (27) together with (39), (60), (61), (62) and (63) gives 



4(c n ^ 2 Cn 4" C n,+lCn c n ,-|-i c n C n C n —\C n C n —\ C n C n —2) 4~ 10LC nJr ^C n 2c n C n _itt — Ijy, 

for n > 2 and 

4(cJ 1 ci + c 2 ~ 1 cic 2 ~ c\ — 6[ cqc± cq) + 2«C2 x ci — 2c{ 1 CQa = Ijy, 
which can be written as 



/3 n+2 /3„+l + /3 2 + i - /3 2 - PnPn-l + Ctfin+l ~ fin® = I 



(64) 



for n > 2 and 



4- Pi - 01 + af3 2 -p ia = I 



N: 



(65) 
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respectively. Using the Abelian character of the set of /3's we arrive to 



Pn+2Pn+l + Pn+l ~ Pn ~ PnPn-1 + "(Ai+l _ Pn) = I 



n = 2, 3, . . . , 



/3 3 /32 + /3 2 2 -/3 1 2 + a(/3 2 -/3i)=I JV . 



Summing up in (66) from i=2 up to i=n we obtain 

n 

X)[A+2A+i + Pl+x - Pi - PA-i + - A)] = X) 1 



V- 



i=2 



i=2 



and consequently we conclude 

fin+2fin+l + Pn+lPn + /?n+l 2 + a Pn+l = ^IjV + &j 

where 

k := /3 2 /?i + /3 3 /3 2 + Pi + aP 2 -l N = P2P1 + Pi + Pi®, 



n > 1, 



where we have used (67). We now proceed to show that k=\jq. (25) implies, for n=l and z = 0, 

^(0) = M 1 (0) Z 1 (0), 



(66) 
(67) 



(68) 

(69) 
(70) 
(71) 



with m[ qS> given in (62). This leads to 



Now, using (|3| we deduce that 



P>(0) 



In 




(0) ( Pi(0) 



c Po(o); • 



Mj 



(o)/0 



Co 



which allows us to immediately deduce that 

P2P1 + Pi + Pxot = h 
and consequently k—Ijsr. Finally, we get 

(3 n+2 P n +l + Pn+lPn + Pn+\ + ^Pn+l = nI N + I N . 

Finally, notice that (74) reads 

P2 = Pi 1 -Pi -a. 



(72) 

(73) 

(74) 
(75) 

(76) 

□ 



This theorem ensures that pi fixes p n for all n > 2, and therefore p n = p n (p\,a). Moreover, 
we will see now that the solutions p n not only commute with each other but also that they can be 
simultaneously conjugated to lower matrices. This result is relevant in our analysis of the confinement 
of singularities. 



15 



3.1 On singularity confinement 



The study of the singularities of the discrete matrix equations (58) reveals, as we will show, that 
they are confined; i.e. the singularities may appear eventually, however they disappear in few steps, 
no more than four. The mentioned singularities in (58) appear when for some n the matrix (3 n is 
not invertible, that is det (5 n = 0, and we can not continue with the sequence. However, for a better 
understanding of this situation in the discrete case we just request that det j3 n is a small quantity so 
that f3~ 1 and {3 n+ i exist, but they are very "large" matrices in some appropriate sense. To be more 
precise we will consider a small parameter e and suppose that in a given step m of the sequence we 
have 

/? m _i = 0(l), det0 m _i = O(l), (77) 

Pm = 0(l), det/3 m = 0(0, (78) 

for e — > and with r < N — 1. In other words, we have the asymptotic expansions 

/V-i = /3 m -i,o + An-i,ie + 0(e 2 ), e -> 0, det /3 m _i, ^ 0, (79) 

/3 m = An,o + /We + O(e 2 ), e^O, dim Ran /3 mj0 = N — r. (80) 

We now proceed with some preliminar material. In particular we show that we can restrict the 
study to the triangular case. 

Proposition 10. Let us suppose that Pi and a are simultaneously triangularizable matrices; i.e., 
there exist an invertible matrix M such that fli = M0iM _1 and a = MjM' 1 with <fi\ and 7 lower 
triangular matrices. Then, the solutions j3 n of (58) can be written as 

f3 n = M<j> n M~ x , n > 0, 

where <p n , n = 0, 1, . . . , are lower triangular matrices satisfying 

4> n+ i = ncj)' 1 - n _i - (f) n - 7. 

Moreover, let us suppose that for some integer m the matrices fim+i, Pm an d ct are simultaneously 
triangularizable, then all the sequence {/3„}^1 is simultaneously triangularizable. 



Proof. In the one hand, from (58) we conclude that M^^M is lower triangular and in fact that 
{M _1 /3 n M} n > is a sequence of lower triangular matrices. In the other hand, if for some integer m 
the matrices /3 m+ i, {3 m and a are simultaneously triangularizable we have 

Pm+l — m fim ~ @m ~ Pm-l ~ a i 

which implies that /3 m _i,/3 m _2 are triangularized by the same transformation that triangularizes 

Pm+i, Pm and a. □ 
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The simultaneous triangularizability can be achieved, for example, when = 0, as in this 

case we can always find an invertible matrix M such that 0i = M$\M~ X and a = M^M~ X where <pi 
and 7 are lower triangular matrices, for example by finding the Jordan form of these two commuting 
matrices. This is precisely the situation for the solutions related with matrix orthogonal polynomials. 



Obviously, this is just a sufficient condition. From now on, and following Proposition [TOj we will 
assume that the simultaneous triangularizability of a and /3\ holds and study the case in where a 
and all the /3's are lower triangular matrices. Thus, we will use the splitting 



n = D n + N n , 
a = dp + cxn, 



(81) 
(82) 



where D n = diag(D n; i, . . . , D n -jsr) and oid = diag(ax> ) i, . . . , old,n) are the diagonal parts of j3 n and a, 



respectively and N n and are the strictly lower parts of f3 n and a, respectively. Then, (58) splits 
into 



D n +i + N n+1 = nD n l - D n _ x -D n -a D 
+nN n - N n _x - N n -a N , 

where N n denotes the strictly lower triangular part i.e., 

= D- 1 + N n . 



(83) 



Hence, (58) decouples into 



N n+1 = nN n - N n _i -N n - a N . 



34) 
35) 



In this context it is easy to realize that there always exists a transformation leading to the situation 
in where 



/ 




An,0 








An,0;r+l,l Pm,0;r+1,2 
(3m,0;r+2,l Pm,0;r+2,2 






Pm,0;r+l,r+l 









m,0;r+2,r+2 



\ 








(86) 



\ Pm,0;N,l fim,0;N,2 " " " Pm,0;N,r+l 0m,O;N,r+2 ' ' ' Pm,0;N,N/ 

Proposition 11. The singularities of the diagonal part are confined. More explicitly, if we assume 
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that (79), (80) and (86) hold true at a given step m then 



D m +2;i 



m -1 a /3m,2;i,i m . ^.i \ 

e — Pm-l,0;i,i ^2 + U \ e )i 



m,l:i,i 



m 



-6 1 + /3 m _l,0;i,i + 



Pm,2;i,i m 



+ 0(6) 



m + 3 2N 



Hi 



(87) 



m/3, 



m—l,0:i,i 



2a D . 



m + 3 



+ 0(6), 



u>/ien e — > 0. 



Proof. From (79), (80) and (86) we deduce 



D m -l,i — Pm-l,0;i,i + Pm-l,l;i,i e 

D m ,i = (3m,i-i,i€ + 0(e 2 ), 



for e — > 0, with i > r + 1. Substitution of these expressions in (84) leads to the stated formulae. For 
i < r the coefficients D TO _i ; j and _D m;i are 0(1) as e — >■ 0, thus they do not vanish, and consequently 
there is confinement of singularities for the diagonal part D n . □ 

In what follows we will consider asymptotic expansions taking values in the set of lower triangular 
matrices 



T:={T + T ie + O(e 2 ), e^O, T* G t N }, t N := {T = (T M ) G C 



WxJV 



Xjj = when i > j}, 



(89) 



where ijv is the set of lower triangular N x N matrices. The reader should notice that this set 
T = tjv[[e]] is a subring of the ring of C iVx7V - valued asymptotic expansions; in fact is a subring with 
identity, the matrix 1 N . We will use the notation 



T; 



T itll 



i > 1, 



(90) 



where Tj n G t r , Tj j22 G and T ij2 i G r ^ xr . We consider two sets of matrices determined by 
(86), namely 



t:={K = ( 
I := {L-x 





-^0,21 -^0,22 

L-x,u 

£-1,21 



^,K ,2lGC^ Xr ,ir ,22GW_ r } ; 



18 



and the related sets 

K := | if = K + Kit + 0(e 2 ) G T, K G «}, 

L := |L = L^e" 1 + L + Lie + 0(e 2 ) G e^T, L_i G [}, 

which fulfill the following important properties. 



(91) 
(92) 



Proposition 12. 1. Both K and eL are subrings of the ring with identity T, however these two 
subrings have no identity. 

2. If an element X G K is an invertible matrix, then X^ 1 G L, and reciprocally if X G L is 
invertible, then X^ 1 G K. 

3. The subrings eh and K are bilateral ideals of T; i.e., L • T C L, T • L C L, T • K C K and 
K ■ T C K. 

4. We have L • K C T. 



Theorem 2. If f3\ and a are simultaneously triangularizable matrices then the singularities of (58) 
are confined. More explicitly, if for a given step m the conditions (79), (80) and (86) are satisfied 
then 



An+l) /3 m +2 G L, 



l3 m +3 e K, 



/3 m +4 G T, 



det/3 m+4 = 0(1), e^O. 



Proof. From (80) and (86) we conclude that (5 m G K and consequently /3" 1 G L. Taking into account 
this fact, (58) implies that (3 m+ i G L. Therefore, G K and (58), as (3 m+ i G L, give /3 m+2 G L 

and consequently /?m+ 2 ^ Iterating (58) we get 

/Ws = f3 m -(m + l)f3~ 1 +1 + (m + 2)/3 m ^ 2 . (93) 

Using the just derived facts, /S^+i; Pm+2 e ^7 an d that An £ K, we deduce /3 m+ 3 G K which implies 
/3 m + 3 G L. Finally, Q gives /3 m+4 as 



/3m+4 = (m + 3)/3 m+3 — /3 m+ 2 — An+3 — 

We conclude that there are only two possibilities: 

1- (3 m +4 = 0(1) for e — > 0, or 

2- Pm+4 ^ L. 

Let us consider both possibilities separately. 



(94) 



1. Recalling that the diagonal part has singularity confinement, see Proposition [TTJ in the first 
case we see that det /3 m+ 4 = 0(1) when e — > 0, as desired. 
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2. In this second case we write An+4 as 
P m+ , = p-l 3 A + 0(l), 



e -»■ 0, 



A := (m + 3)1 - (3 m+3 (3 m+2 . 



Observe that repeated use of (58) leads to the following expressions: 

A = I+[(m+ l)(5 m l +1 - An] An+2 

= k + I-[{m+ - An]An+l 

= k - ml + (3 m (3 m+1 

= k - An(An + fim-i + a), 



where 



k := \{m + - p m ][(m + l)/3j +1 - An - a]. 



(95) 



From these formulae, as An +1)/ S m e ^ we deduce that k e K and also that (3 m ((3 m + fi m ~\ + 
a) G K. Hence, we conclude that A G K and from (95) and Proposition 12 we deduce that 
An+4 = 0(1) when e — > 0. Consequently, we arrive to a contradiction, and only possibility 1) 
remains. 

□ 
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